Spherical Universes with Anisotropic Pressure 



J R Gair 

Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CBS 
OHA, UK 

E-mail: jgair@ast.cam.ac.uk 

Abstract. Einstein's equations are solved for spherically symmetric universes 
composed of dust with tangential pressure provided by angular momentum, L{R), 
which differs from shell to shell. The metric is given in terms of the shell label, R, and 
the proper time, t, experienced by the dust particles. The general solution contains 
four arbitrary functions of i? - M{H), L{R), E{R) and to{R). The solution is described 
by quadratures, which are in general elliptic integrals. It provides a generalization of 
the Lemaitre-Tolman-Bondi solution. We present a discussion of the types of solution, 
and some examples. The relationship to Einstein clusters and the significance for 
gravitational collapse is also discussed. 



PACS numbers: 0420-q,0420Jb,9880Hw 



Submitted to: Class. Quantum Grav. 



Spherical Universes with Anisotropic Pressure 



2 



1. Introduction 

In 1939, Einstein solved the field equations of general relativity to describe a cluster 
of particles moving in circular orbits about the centre of the system. The static and 
spherically symmetric solution was supported by a balance between gravity and the 
centrifugal force due to the angular momentum of the particles. These Einstein clusters 
have been studied extensively. Gilbert |^ investigated their stabihty, Florides [§] took 
one as a model for the interior of a star and Comer and Katz 0] used them as a model 
of a 'heat bath' for a black hole. A natural generalization is to consider a similar 
system which evolves dynamically. This was first considered by Datta [^], who limited 
his discussion to the case in which all the particles have the same angular momentum. 
Bondi generalized Datta's model to allow each shell of particles to have a different 
angular momentum, and discussed the types of behaviour that might occur. The same 
model was discussed by Evans 0, who investigated the system for both thin heavy shells 
of particles, and a cloud of dust. None of these authors obtained a general solution for 
the metric of the spacetime. Datta obtained the solution for the simple case in which the 
angular momentum is the same for each shell. Bondi considered matching the general 
solution to an external Schwarzschild metric, but did not find any additional solutions. 
Evans discussed how the system would evolve from initial data, but again did not obtain 
any solutions. 

In recent years, much work has been done on gravitational collapse, as an 
investigation of the Cosmic Censorship Hypothesis (see for a recent review and 
references therein). In this context, Magli examined this model again. He discussed 
a general class of spherically symmetric solutions to Einstein's equations with vanishing 
radial pressure. He was able to solve for the metric in terms of mass-area coordinates 
10| , based on the method used for charged dust collapse by Ori [|lT|. Using his solution, 



Harada et al [|12| investigated naked singularity formation in the model, and gave a 
particular case in which the metric could be written in terms of elementary functions. 
This special case was investigated by Harada et al and Kudoh and singularity 



formation in the general case was examined further by Jhingan and Magli |T5[ . 

We tackle this problem in a new way, solving for the metric in terms of the time 
experienced by the dust particles composing the universe. The mass-area approach has 
the advantage that the solution reduces to a single integral which is in terms of the 
coordinates used. This is not true in general in the new approach. However, the mass- 
area solution makes the structure of spacetime unclear, since both these coordinates can 
be spacelike in general. Moreover, such coordinates are only appropriate to monotonic 
collapse or expansion, in which the area can be used as a time coordinate. The mass- 
proper time solution is better suited to studies of model universes, which may have 
phases of both expansion and collapse, and it is physical, since the time coordinate is 
understandable, as the time experienced by the dust. 

In section 2, we describe the model and derive the equations. In section 3, we solve 
the equations, present a general discussion of the types of solution, and describe the 
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relation to Einstein clusters. In section 4, we present some examples. Section 5 gives 
our conclusions. 

2. The Model 

The spacetime is composed of a large number of dust particles, which can move radially, 
and have angular momentum about the origin. Spherical symmetry is assumed, so that 
all quantities are functions only of the time and the radial coordinate. To produce this 
symmetry, at each point there must be a large number of particles, moving with the same 
angular momentum, but in every tangential direction. A comoving radial coordinate, 
R, is used, i.e. particles remain on surfaces R = constant for all time. The general 
spherically symmetric metric is given by the ansatz: 

ds2 = e^-^^^'^^di^ - e2^(^'*)di?2 - r^R, t){de^ + sin^ M^^). (1) 

We assume that the dust only interacts with itself gravitationally, so that each 
dust particle moves on a geodesic of the spacetime. As a consequence of the geodesic 

equations, the angular momentum, L — r'^\J (d6'/dr)^ + sin^ 9 (d0/dr)^ is a constant for 
each particle (where r denotes the proper time for the particle). Spherical symmetry 
requires that L be the same for each particle on a given shell, i.e. L — L(R) only. 
The energy momentum tensor for dust is T"'^ = pu"u^, where p is the proper density 
of the dust and = dx'^/dr is the 4- velocity. This would be T'^^ for this spacetime 
if the 4-velocity of the dust was well-defined. However, particles are moving in all 
tangential directions at each point, and so the 4-velocity is not single valued. But, as 
the particles are assumed not to interact, the overall energy-momentum tensor can be 
taken as a sum of the individual contributions over all the particles T'^^ = ^ pu^u". 
The choice of comoving coordinates makes the radial pressure —T^ vanish (u^—O). 
Spherical symmetry of the overall system requires the two tangential pressures to be 
equal Tq — — —p{R,t). For a single dust particle, UfU^ — e^^ldt/dT^. Prom 
the geodesic equations, c^'^ (dt/dr)^ = 1 -|- L'^r"^. Since L is conserved, we have 
u^ug + u'^u^ = —L'^/r^. The ratio of u^ug + u'^u^ to UtU^ is the same for every particle 
on a shell, and so the energy momentum tensor is given by 

MR, t) = -Tg' - = ^^T/ = j^AR^ (2) 

e{R,t) is the energy density, T^, which includes a contribution from the kinetic 
energy of the particles. 

Conservation of the energy- momentum tensor, T^.^, yields the two non-trivial 
equations: 

r' 

i/e - 2-p = (3) 
r 

e + Ae + 2^(e-Fp) = 0. (4) 
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In these, ' = {d/dR)^ and ' = {d/dt)^. Using the convention -R^^^ = T^^^ — 
Qap _ j^a/B _ Ig^f^ji for the Einstein curvature tensor and the identification r = e^, the 
non-trivial Einstein equations are given by: 



-TA^2e- + A,' -.',)^0 (5) 

^T/ = e-- (2V + if^r) + e-''^ - e^^ {2^^" + 3 i^^'f - 2AV') - A = (6) 
^""^Tl = e-^'^ (2/. + 3 (/i)' - 2z>/:i) + e'^^ - e^^^ ((/x')' + 2z/V) - A = (7) 



c4 



SttG 



Tg^ = e ^Az> + z>/i — A/i — A — A^ — /i — /i^ 



I -2A ( II , /2 I // I /2 \/ / /\/ I / I A SttG . . 

+e (z/ +1/ +/i -A;U-z/A+z//ij+A = -^P- (8) 



These are the field equations of Landau and Lifshitz ([|T^ §100), but we have written 
2A for their A, 2yU for their /i and 2v for their v and included the cosmological constant 
A. In the rest of the paper the areal radius r = will be used. From (H) and (^, 

1^2 ^1 

Substitution into equation (||) and integration gives 



2A ^ V .^Q^ 

1 + 2E ^ ^ 



In (10), the constant of integration has been taken to be — ln(l + 2E{R)). In the 
Tolman-Bondi case, it is usually denoted as 1 + f{R)- However, the function 1 + 2E{R) 
should be regarded as the relativistic energy, i?^, and in the Newtonian limit becomes 
1 + 2E-^{R), as will become clear. The use of E{R) instead of / makes this identification 
more explicit. Substitution of and (|r^) into (|^ and integration gives: 

1 _ GAm _ 1^^2 ^ ^ _Er^ 

2 r 6 2(L2 + r2) + r"^' ^ ' 

The constant of integration has been taken to be —2GM{R). On substitution of 
(I), (0) and dn]) into (D we obtain 

2GM'(R) 8nG /"^ . 9 e , , ^ 

^ ' - ' - M{R) = / 47rr^—r'dR. (12) 



r^r' Jq c^ 



M{R) is the Misner-Sharp mass fl^, which is conserved for spherical systems with 
vanishing radial pressure. This is not the energy density summed over all shells. For 
that, the factor exp{X) must be included in the integral to give the correct volume 
element. 
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3. Mass-Proper Time Solution 

The equations (|^) and (|Tl|) are coupled equations for the two remaining unknown 
functions, r{R,t) and v{R,t). They cannot be solved explicitly in the coordinates t 
and R. Magli found a solution by using mass-area coordinates, that is r and R. 
We shall obtain a solution in terms of a different time coordinate. There are two times 
besides t relevant to this problem. The time experienced by an observer at rest on a shell 
obeys dT^ = e^'^dt^. The equation of motion in T is therefore identical to ([TTl) , with the 
first term replaced by \{dr/dT)\. Alternatively, there is the proper time experienced 
by a dust particle, r, which obeys 



On changing coordinates from (t, -R) to {t,R), equation ( ]TT| ) becomes: 

1 fdrV GM f 1 „ 1 „ 1L2 

: =E + 1 + — +-AL2 + -Ar2-- — . (14) 



2 \dr J ^ r \ r"^ J 6 6 2 r"^ 

This equation is exactly the equation of motion of a test particle with angular 
momentum L and energy + 2E in a vacuum Schwarzschild metric of mass M. This 
result is a consequence of the assumption that the dust particles follow geodesic paths, 
and justifies the identification of E as the energy and M as the mass made earlier. 

Equation (ITl) may be integrated, keeping R constant: 



dr 



To{R) is a constant of integration. The function A{r, R) is given explicitly by ([I^) , 
since dr/dr is known as a function of r and R. In the case A = 0, the integral is an 
elliptic integral. For non zero A the integral is hyperelliptic. A{r, R) = t — tq{R) gives 
the function r{R,T) on inversion. We obtain [dr/dR)^ from ([T5|) as 

The solution is not yet complete, as the metric has not been written in the new 
coordinates. To do this, dr must be expressed in terms of dR and dt. Equation (0) 
gives (9r/(9t)ij. To obtain {dT/dR)t, we first differentiate (^3|) with respect to R and 
then change the order of partial derivatives, obtaining an expression for {d / dt){dT / dR)t- 
This may be written as an equation for {d / dT){dT / dR)t by noting 
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The resulting equation has an integrating factor 1 + L r , giving the result: 

m)^^~L^ + r^ {m)^/ r2(9r/9r)^^'') ■ ^^^^ 

Again, this is an elliptic integral if A = and hyperelliptic otherwise. It is 
closely related to the integral for r{R,T). Evaluation of the integral, (19) includes a 



constant of integration, go{R) say. This defines the origin of proper time, which may 
be different on different shells. This function is eliminated by choosing a common 
origin of time for all shells, i.e. at some time, say t = 0, we set r = for all R. 
Knowing the distribution of matter at that time, r(0,i?), the equation {dr / dR)^^^ = 
with r = r{0,R) determines go{R). This does not remove any generality from the 
problem, as the proper time is only defined locally for a shell. Adjusting where the time 
is measured from does not physically affect the spacetime. Using (pIS]), ([19|) and the 
identity e^'dt = \/l + L'^r~'^ (dr — (dr/dR)tdR), the metric, (|l|) can be written 

= f 1 + dr^ - 2 f 1 + ^ drdi? - r\d9^ + sin' 9d^') 



r"^ J \ r"^ J \dR J ^ 

r' + L'^ ( f dr\^ ( d 



r 



2( 




^1 +21^ 



1 + 2E) \ \dRj^ \dR 




All the functions in (|20|) are now known. r{R, r) is given by inversion of (|15|), {dr/dR)^ 
is given by (|16[) and {dr/dr)^ by (|T^). The function (dr/dR)^ is given by (|1^), noting 

dBir,R) \ _(dB\ ^(dr\ ( dB_\ 
The metric (pOD , with the equation of motion ([l^) , the Misner-Sharp mass ([l2|) and 



expression (19) for {dT/dR)^ solves the problem in terms of the proper time and shell 
coordinate. The evolution of the system depends on the four free functions M, E, L 
and tq. One of these corresponds to a choice of the shell label, R. The others are fixed 
by the initial density, kinetic energy and angular momentum of the dust. 

3.1. Shell Evolution 

Equation ([T^ ) determines the motion of each individual shell. It may be written more 
clearly as: 

Equation (|22|) is the familiar Newtonian equation except for the factor (1 + L^/r^). 



This is a relativistic correction, which is due to the effective mass of the angular 
momentum. A physical evolution must have E — V > 0. The shape of the potential V 
determines which types of motion are possible, and the energy determines which regions 
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of the potential a shell may penetrate. In the case A = 0, the potential was discussed 
by Bondi 0. 

Figure la illustrates the possible shapes of V"(r, i?) if A = 0, which is familiar as 
the potential for a Schwarzschild solution of mass M, angular momentum L and energy 
Vl + 2i?. If A is non-zero, V contains the additional term — |Ar^ (1 + L^/r^). This 
has little effect at small r (only a constant shift — |AL^), but dominates at large r. 
For positive A, it causes the potential to turn down towards — oo, which can add an 
extra turning point in r > 0. For negative A, it causes the potential to turn up, and 
V +00. A single example of each sign is shown in figure lb. The horizontal lines in 
figure lb represent different choices of the energy E. The physical regions for a given 
value of E have V < E, indicated by the solid portions of the lines. The evolution of a 
shell depends on whether the energy intersects the potential above or below the turning 
points. 

There are six possible types of evolution, which are labelled in figure lb. 

(i) Expansion and recollapse. 

This occurs if ii^ — ^ > for < r < Rmax with V{Rmax) = E. A shell in this 
region expands from a Big Bang at r = 0, to its maximum radius, Rmax, and then 
recollapses. 

(ii) Bouncing Universe 

This occurs if V{Rmin) = E and E — V > for Rmin < r < oo. This represents a 
shell which collapses from some initial size, reaches a minimum radius, Rmin, and 
then expands to infinity. 

(iii) Expansion to/collapse from infinity. 

This occurs if E — V > everywhere. A shell starting at infinity can collapse to 
the origin r = 0, or a shell created at r = can expand to infinity. 

(iv) Oscillating Universe 

This occurs if V{Rmin) = E = V{Rmax) with £" - 1^ > for Rmin < r < Rmax- A 

shell in this region oscillates between the two extremes, Rmin and Rmax- 

(v) Circular Orbits 

These occur when a shell sits at a minima (v.a) or maxima (v.b) of the potential. 
Orbits at maxima are unstable, those at minima stable. 

(vi) 'Coasting' Universe 

A shell will coast if the potential has a repeated root at Rare, but the shell is 
initially at r 7^ Rdrc- A shell starting at i? > Rdrc will collapse, but the rate of 
collapse decreases as it gets closer to Rdrc, so that it does not reach it in a finite 
time. A shell starting at i? < Rdrc and initially expanding (perhaps from a big 
bang at r = 0) will coast outwards to the limiting value. The limits in each case 
are the circular orbits of (v). 

If the energy in (vi) is adjusted to E + 6E, the result is a 'hesitating' universe. 
This is one which expands, slows down ('hesitates') as it nears the turning point, before 
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either undergoing accelerated expansion again {5E > 0) or recollapsing {5E < 0) (see 
p!8| for a description of this in the FRW context). 




Figure la Possible potential shapes for A = 0. As the parameter k = L/GM changes, the potential changes 
(starting with the lowest, solid curve) from having no roots and no turning points {k < 2\/3); no roots and a point of 
inflection (k = 2-\/3); two turning points and no roots (2-\/3 < fe < 4); one repeated root (k = 4) and finally two roots 
(fc > 4). 




20 



25 



30 



35 



Figure lb Shape of potential for A > (soUd line), A = (dashed line) and A < (dot-dash Une). The 
horizontal lines correspond to different choices of the energy E. Physical regions have E > V, and are indicated by the 
solid parts of these lines for the A = potential. 
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3.2. Global Behaviour 

The universe is composed of an infinite number of shells of the type discussed in the 
previous section. These shells lie within one another, and so care has to be taken when 
constructing a global solution. In general, a physically reasonable solution is one which 
starts from reasonable initial data - either a Big Bang type singularity or a regular 
matter distribution, and evolves without the formation of shell crossing singularities. 



Behaviour of the Functions of R There are four free functions of R in the solution. 
The time function tq{R) specifies the initial distribution of matter. In big bang models, 
r = for all R at t = 0, which we take to be r = on every shell. In collapse models, it 
is sensible to choose r = i? at t = 0. The mass function M obeys equation (0). If the 
spacetime is to satisfy the energy conditions, we require e > and so M' and r' must 
have the same sign. Finally, in studies of gravitational collapse, it is conventional to 
require the initial data to be regular and all functions to be Coo until a singularity forms. 
This leads to the conditions M{R) = M^R^ + M^R^ + E{R) = E2R^ + E^R"^ + ... and 
L'^{R) = LiR'^ + Lg-R® + ... (see |jT^). In the other models, for instance the expansion 
and recoUapse case, the condition of C^o evolution is not necessarily relevant, as the 
universe is born from a singularity. 



Curvature Singularities A curvature singularity occurs when an invariant diverges, 
for instance R'^. Expressions (D and give R'^ = - (87rGr°/c^ + 4A) = 

— (167rG'^M'/(c'^(L^ + r^)r') + 4A). This is singular if r'/M' = 0, which corresponds 
to shell crossing, or when r = provided L = there. In ||12[, they took L oc R^, in 
which case the centre can become singular when the central shell R = reaches it. As 
further shells reach r = 0, the point remains singular since there is a non-zero point 
mass there. In general, the central singularity can be naked or clothed (i.e. a black 
hole). It can be naked only for the time between the shell R = reaching r = and any 
subsequent shell arriving. Photons cannot escape from any region with r < 2CM{R). 
For the central shell, which has M = 0, photons may escape from the origin, but this 
will not be the case for shells with M ^ 0. This was also mentioned by Harada et 
al. The formation of naked or clothed singularities will be discussed for some specific 
examples later. 



Coordinate Singularity at E = —1/2 It is clear from the metric (20) that if E = — |, 
the guR coefficient is singular, unless r' is also zero at that point. This is a coordinate 
singularity and may be removed by a coordinate transformation R —>■ R' with dR' /dR = 
1/VI + 2E near the singularity. In the new coordinates, {dr/dR')t = at this point, 
but this does not represent shell crossing as dM/ dR' = also. This is the equator of the 
universe. On the other side of the singularity, the universe is in reverse, with M' < 
and r' < 0. The other half of the universe may be a mirror image of the first, or of 
any other model which has the same values of M, L and E at the equator. In this way. 
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two half universes are joined together at the coordinate singularity. This is similar to 
the k = —1 FRW model (^). The point R = 1 is singular, but writing R = sinx, this 
singularity is removed. The universe consists of two identical hemispheres, joined at the 
coordinate singularity R = 1. 



Shell Crossing The solution breaks down when a shell crossing takes place, which is 
indicated by r' = and the density, e, diverging. A similar problem occurs in the 
Tolman-Bondi solution (with L = 0), and has been discussed in detail in |jl9l. In more 



general solutions of relativistic stellar dynamics (e.g. [^), such problems do not 

arise. These models use the Einstein- Vlasov equations, the evolution is followed in 
phase space and streams may cross. However, the non-static Einstein- Vlasov system is 
difficult to solve. In fact, this model could be regarded as Einstein- Vlasov, because the 
dust is assumed to follow geodesic paths. The phase space density is proportional to 
6{pr — r{r,t)) at each point. There is no radial velocity dispersion and hence no radial 
pressure, but the phase space density is divergent, since the momentum is confined to 
lie in a plane at each point. In more general Einstein- Vlasov systems, there will be some 
radial velocity dispersion. This gives rise to radial pressure, which has been ignored in 
this model and is likely to be important at shell crossing. In fact, the problem at shell 
crossing is not with gravity, but with the solution. The shell crossing problem could be 
avoided by looking at the evolution in phase space, and indeed the divergent density is 
only a surface layer, which can be treated in general relativity (|2^). Our method of 
solution, however, relies on the fact that the shells are distinct, which gives an equation 
of motion for each shell that is independent of the others. Moreover, when the density 
becomes large, interactions between the dust particles are likely to become important. 
Thus, this model shall only be regarded as valid until a shell crossing first takes place. 

A shell crossing is indicated by {dr/dR)t = 0. This derivative is taken at constant t 
rather than by holding one of the other times, r or T, constant. That this is the correct 
derivative is clear from ([l^), which tells us that e oc {dr / dR)^^ . When r' = 0, the density 
diverges as the shells attempt to occupy the same space. This is not true if M' = r' = 0, 
which may occur at the equator of a closed universe, as discussed previously, or where 
e = 0. The latter indicates a region in which the shells are unoccupied and hence shell 
crossing is not a problem. With this exception, the condition for no shell crossing is 
that 

(^) +(^) (^) ^0 V r,R. (23) 

R 



dR , \dR \dR ,\dT 



The functions {dr/dr)^ and {dT/dR)^ are given by ([1^) and (0) respectively. 
(dr/dR)^ may be obtained from ( pT]) or by differentiating (p!^ and then evaluating the 
expression: 

fdr\ r-, r f GM'r^ + (GML^)' - LL'r + E'r^ + \ALL'r^ , , 
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Certain choices for M, E, L and tq will give solutions in which shell crossing never 
occurs. Some necessary conditions for such an evolution will now be discussed. In the 
following, M' > is assumed, and so physically r' > is required. Regions with both 
of these negative are equivalent to this case by relabeling the shell coordinate R. 

If a solution begins with a big bang at t = 0, then equation (0) tells us 
that for r <C ri (where ri{R) is the location of the first turning point of V{r,R)), 
r ~ (5r)^/^(G'ML^/2)^/^. Integration of ([T9| ) with the boundary condition t = = t 
gives {dT/dR)^ = -SL/i^ (L/zi^)' r/(L2 + /iVt) where = (25GMLV2)1/io. Using 
these in (p3D, the condition for evolution without shell crossing becomes (In {fi^/L))' > 0, 
but this just states M' > 0, which is assumed. 

In unbounded evolution {dr/dT)'^j^ ~ |Ar^ + 2E + |AL^ asymptotically. If A 7^ 0, 
this is the same for all shells. If A = equation (p3D becomes {dr/dR)^ ~ A + for 
large r. If E' < 0, shell crossing will inevitably occur. We therefore require E' > for 
unbounded expanding evolutions. For a collapsing evolution that begins at infinity, the 
sign of the E'/E term changes and the reverse must be true, E' < 0. 

If any shell has a maximum radius, all the shells within it must. Similarly, if a 
shell has a minimum radius, all those exterior to it must. Moreover, the maximum or 
minimum radius must be an increasing function of R. This applies to the maximum 
radius in the expansion/recollapse case, the bouncing point for a bouncing region and 
the limiting radius in a coasting evolution. 

In a bouncing region of the spacetime, the bounce point, i.e. the part of the 
universe that is bouncing, must move outwards 'supersonically', that is faster than the 
surrounding shells expand. In a dust collapse model, it has been shown for the Lemaitre- 
Tolman-Bondi case by Hellaby and Lake among others that it is always possible to 
choose the initial conditions in such a way as to ensure there is no shell crossing. The 
freedom in tq allows the same here, although equally initial conditions can be chosen 
which ensure shell crossing does occur. 

Finally, for a coasting region of the universe (case (vi)), {dr/dr^j^ ~ f{R){rQ{R)—r) 
for every shell asymptotically. Equations ([liD , (|T^) and (|23|) tell us that (dr/dR)^ ~ 
r'o + T{f'{R)-L'^f{R){ln{L/r^)y/{L^ + r^))exp (-/(/?)r). If dro/di? and the bracketed 
term are both positive, there will be no shell crossing asymptotically. This does not 
exclude shell crossing at earlier times, however. 

A family of self-similar solutions to this model can be found. In that case, the shell 
crossing conditions can be written as inequalities for the free constants in the solution. 
The self-similar models will be discussed in a separate paper. All the types of evolution 
(i)-(vi) can occur for the whole universe without shell crossing in the self-similar case, 
except for (iv), the oscillating universe. 

In the oscillating case, shells will be continually changing directions for all time 
and so shell crossing seems quite likely. In the next section a solution undergoing small 
oscillations will be derived by perturbing an Einstein Cluster. If a solution undergoing 
large oscillations is to exist, a necessary condition is that the period of oscillations in t 
is the same for all shells. This is difficult to evaluate in general, as an explicit expression 
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for is not known. Progress can be made using (|23|). If the system is undergoing 
oscillations, then necessarily r(i?, T + nI\.{R)) = r{R, T),'^n G Z, i.e. each shell oscillates 
with some period n(i?) in r. If we denote the limiting radii of the oscillations of each 
shell by rj„(i?) and rout{R), then this period is given by the integral: 



U{R) 



rout{R) 



-dr. 



(25) 



irUR) VAr^ + (6E + AL2)r4 + QGMr^ - SL'^r'^ + GGML^r 
From (p3|) and (p!9D , it is possible to relate the shell crossing condition after n 
periods have passed to the values of the functions in the first period: 

Lr^ f dr\ f d 
dR. 



--To+nn{R) 



dr \ 



dR L'^ + r^ \ dr 

/to V 



T"0 



+ n 



dr 



TO 



L2 



-n'(i?) -2L 




(26) 



This uses the fact that {d/dR)ro+nn{R) = {d/dR)ro - nU'{R){d/dT)ro. Taking 
< To < n(i?), the first line of (^61) is the value of {dr/dR)^ on the first cycle. The 
second line is proportional to {dr/dT)^^^^, which changes sign during a single oscillation 
and n, which may be arbitrarily large. Thus, even if shell crossing does not occur during 
the first period, it will occur inevitably for large enough n, unless the term in square 
brackets is 0. This is a function of R only and gives the condition: 

rout(R) T 



d 



r-out(-R) ]^ 

rUR) ^2{E-V{r,R^ 



dr 



''^2{E-V{r,R)) 



dr 



0. (27) 



This is the statement that the period in t of every shell is the same. It is a necessary 
condition on the functions M, L and E for the evolution to proceed without shell 
crossing. The condition is only sufficient if there is no crossing during the first oscillation. 
The additional freedom in the initial data allows the construction of oscillating universes 
which have no shell crossing on the first cycle and hence never. 



3.3. Einstein Clusters 

If every shell is in a circular orbit the spacetime is static. In a circular orbit f = f 
i.e. V{r, R) = E{R) and dV/dr = 0. Labeling the shells by i? = r this gives 



0, 



L2 



2E 



-1 + 



IAR2 _ 2^)2 



l~3^ ' 1-3^ 



A circular orbit is stable if 

'd'^V\ M(l-6M)+Ai?2(6M_| 

\r=R 



> 0. 



(2^ 



(29) 



dr^J,_„ i?2(l-3M) 
For this static case, the integrand in ( p^D is independent of r and so {dT/dR)t is 



given easily. Alternatively, equation 

M 



V 



\KR- 



R{\ - 2M - |Ai?2) 



gives 
di?. 



(30) 
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The solution in (r, R) coordinates is given explicitly by specifying M and A and 
in (t, R) coordinates by specifying v and A. However, quadratures are required to 
determine M oi v from the matter distribution e. 

Zero Cosmological Constant For A = 0, this solution is the Einstein Cluster 
Expression (|30D and the stability criterion (|^), < M < |, agree with those derived 



in a more roundabout way in 0. Expression ( pSD requires M < | if L is to be real. 
The case M = | corresponds to a universe composed of radiation rather than dust and 
will be discussed elsewhere. This spacetime has —j^<E, which means there is never 
enough mass to close the universe (i.e. E can never equal —1/2). 

Non Zero Cosmological Constant If A > 0, but L = 0, the solution is the Einstein 
Static Universe, which has M = |Ai?^ everywhere. The energy E = — |M = — | when 
M = |, so the universe is closed with a maximum radius Rmax = 1/V^- From 
this universe is unstable everywhere. If both A and L are non zero, changes sign 
when M = I and when M = |Ai?^. One of these equations always gives a maximum 
radius for the cluster. This contrasts with the A = case in which infinite clusters are 
possible. If M = I = |A_R^ at the same point, L'^/R^ will be finite there, and the energy 
E reaches — |, the energy required to close the universe. The static model can therefore 
represent either a cluster or a closed static universe. In the latter, the maximum radius 
of the universe, Rmax = is independent of the L distribution. 

From ( p9D , the universe is stable if x — ^ ^/y < M < x + ^/y, where x = + 1 Ai?^ 
and ?/ = 1 - 28Ai?2 + SGA'^R^. This is satisfied near R = if |Ai?2 < but it never 
holds for M > |. The outer regions of a closed static universe must always be unstable, 
which is the case everywhere in the Einstein Static Universe. 

The A < case is very similar to A = 0. Only cluster models are possible, 
which may be stable or unstable, and M < 1/3 so that L is real. The stability 
criterion is slightly relaxed. The cluster is still stable where < M < 1/6 and 
is unstable for M > 5/18. For 1/6 < M < 5/18, the cluster is stable if i? > 

(M(6M - l))/(|A|(5/3 - 6M)) and unstable otherwise. 

Cosmic Censorship The Einstein clusters are the asymptotic states for 'coasting' 
evolution. Such evolutions have the same E, M and L, but the evolution begins away 
from the circular orbit. In the A = case, marginally bound coasting collapse was 



discussed in |12]. This solution has E = and hence L = AGM (from (|28|) ). The 
solution exhibits a globally naked singularity at the origin, which is partially due to 
the presence of a naked singularity in the asymptotic Einstein cluster. For an Einstein 
cluster with a singularity at the centre e ~ R~^ and M oc R^'^ near R = with 
< 5 < 2 (so that M is bounded). The expressions (0) and (|30|) allow the radial 



photon equation dR/dt = exp(z/ — A) to be evaluated. For < 5 < 2, we find 
exp(z/ - A) ~ (1 - 2MoR'^-^Y^-^^^^ near R = 0, where Mq = lim^^o(M/^^"'^) is a 
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constant. In the case 5 = 2 this becomes exp(z/ — A) ~ ^2Mo/{i-2Mo)_ either case, the 
photon equation is integrable, photons can propagate from the origin in finite time and 
the singularity is naked. This excludes the marginal case 5 = 2 and Mq = |. This is 
the radiation limit, and will be discussed elsewhere. 

All these Einstein clusters exhibit locally naked singularities, which may be globally 
naked depending on the structure of the rest of the spacetime. However, if the Einstein 
cluster is to be formed as the endstate of coasting collapse, the orbit of each shell in the 
spacetime must be located at a maximum of the potential, i.e. at an unstable point, 
so that E — V > Q in a. region about that point. If A = 0, this requires M > i. This 
rules out all the cases with < 5 < 2, which have M — > as i? ^ 0. Only if 5 = 2 and 
Mo > I can the cluster form in collapse. This leaves a one parameter family of collapse 
solutions which do have naked singularities in the endstate, characterized by the value 



of M in that endstate. The case considered by Harada et al [0 is one such example, 
which has Mq = \- The presence of a non-zero A does not affect the solution near 
r = 0. Again, a locally naked singularity will be present if M ~ Mq > I, a constant, 
near i? = 0. This may be globally naked, depending on the structure of the spacetime 
away from R = 

It is possible in this way to produce naked singularities readily in this model. 
However, fine tuning of the functions L and M is required, and if the solutions form 
in collapse, they must be at maxima of the potential. For that reason, the solutions 
are unstable, since it is likely that even spherically symmetric perturbations will cause 
the cluster to evolve away from the static solution, resulting in rapid clothing of the 
singularity. 

Oscillating Einstein Cluster Perturbation of a stable Einstein cluster produces a 
solution undergoing oscillations. Keeping M and L unchanged, but letting r — > 
R + 5r{R, t), pq{R) + 5p{R, t) and E Eo{R) + 5E{R), we obtain from (|lj 



d^5r 1 9V M(l - 6M) + AR\6M 



5^ 



dr^ 2 dr^ \r=R R2{1 - 3M) 

{d'^V/dr'^) > necessarily, as the cluster being perturbed is stable. This 
is the simple harmonic motion equation - each shell undergoes SHM with period 



CO = -^(M(l -6M) + Ai?2(6M- |))/(/22(i _3M)) in r. To avoid shell crossing, 
the period in t must be the same for each shell. To lowest order, the period in t is 
Q = a;e'^''(l + L'^/R'^)~^. Requiring this to be constant we have 

M' . . '-^ = -. 32 

^ 2M(l-6M + |Ai?2)(i _3M) R 



In the case A = 0, this is integrable and R = (VM(1 - 2Mf''^)/{yt{l - 'iMf'^). 
This reaches a maximum at M^ax = (1 ~ a/2/3)/2, which means there are two solutions, 
one with < M < M^ax and one with M-^ax < M < 1/6. The functions R{M) and 
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M{R) for the two solutions are shown in figure 2. For the second solution, dM/dR = 
before Rmax is reached. The solution must be cut off there, to ensure e > 0. It can be 
matched smoothly onto an exterior, since the density reaches at the boundary. 

This solution is approximate, since terms of order 5^ and higher have been ignored. 
These terms will be important at late times, which will lead to shell crossing eventually. 
The full condition (P7D can be used to compute the SE{R) necessary to avoid shell 
crossing completely. 




(a) G.../R (b) 

Figure 2 (a) Plot of R{M) for the A = vibrating Einstein Cluster. The curve has a maximum at 
M = Mmax = 5(1- ■ (b) M(R) for the two branches of (a). The solid line is the branch with < M < Mmax 

and the dashed line is the branch with M^ax < M < i. 



3.4. Matching 

This model may be used to represent a region within a background spacetime. In this 
context, it is useful to know how to match this solution onto other metrics. 



FEW If we take M = po^^ L = 0, 2E = kR^ and r = Ra{t), then this solution 
is the FRW metric with scale factor a{t): 

ds^ = dt^ - a" it) ( Wdi?2 ^ ^2. ^^2 ^ gi^2 ^^^2^^ ^ (ggN 

\l + kR^ J 

This may be open (A; > 0), fiat {k = 0) or closed (A; < 0). If the functions M, E and 
L are chosen to have this form asymptotically, the solution matches smoothly onto an 
FRW exterior. Such a model represents an FRW universe containing an anisotropic 
pocket which may recoUapse to form a black hole. Such a 'swiss-cheese' universe, 
having Tolman-Bondi regions within an FRW exterior has been considered previously. 



for instance in |]23[. In this model, an example of a one-cell universe is to choose 
2GM = R^, 2E = kR^ and L = R^{R - if for i? < 1 and L = for R > 1. This 
solution does not exhibit shell crossing for k = ±0.25 and A = —1.0, 0.0 or 1.0. In 
each case, the universe is born in a big bang and the core recollapses to form a black 
hole. The evolution of the rest of the universe is determined by k and A, and it either 
expands forever, or recollapses in a big crunch. 
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Schwarzschild-de Sitter If the matter distribution is finite in extent, then spherical 
symmetry requires the exterior vacuum metric to be Schwarzschild-de Sitter: 

dT^ 

ds^ = F(r)dt^ - - r\de^ + sin^ M^'). (34) 
F[r) 

In this, F{r) = 1 — 2Gm/r — |Ar^. The matching of this solution onto a 
Schwarzschild-de Sitter exterior was discussed partially in . We use the Israel-Darmois 



2^ matching conditions for the match. If the density is zero at the edge of the interior 
region, M' = = L at the boundary. The match can be made smoothly by taking 
M = constant and L = for the whole exterior region. If the density is non-zero at 
the boundary, there is only a jump discontinuity in the energy tensor there. This 
is classified as a boundary surface, over which the first {gij, the induced metric) and 
second {Kij, the extrinsic curvature) fundamental forms are continuous. We denote the 
interior coordinates as R,T,r{R,T),9,(f) and the exterior coordinates by t,r, 9,(j) and 
take the boundary to be at i? = Ro, r = ro{t). Comparing the metrics induced on the 



boundary by the interior solution (20) and the exterior solution (0), and choosing the 
angular coordinates 6 and to be continuous, continuity of the first fundamental form 
requires that: 

fo{t)=r{Ro,T) (35) 



dr 
dt 



1 _ 2GA/ _ 1 A^2 

^ — (36) 



L2 



In (|3q) , r = fo{t) = r(r, i?o). Continuity of the second fundamental form 
ensures that because the particles in the shell undergo geodesic motion of the interior 
metric, their motion will be geodesic of the exterior as well. Geodesies of the 
metric (|34D are characterized by two conserved quantities, the energy E and angular 
momentum L, which are the same for every particle on the boundary by spherical 
symmetry. The geodesic equations then give an expression for the evolution of the 
boundary. A second expression is derived from the interior solution, since dr(r, Ro)/dt = 
{dr / dt){dr / dr) and these are given by (|36|) and (|^). Equating the two expressions 
for (dfo/dt)^/F^(fo), we find 

fl) ffi (M(B«) - m) (l + a?^) + 1 + 2B(B„) 

It is clear that the two sides of equation (|37|) are consistent only if L = L{Ro), 
m = M{Rq) and E"^ = 1 + 2E{Rq). With this identification, the evolution of the 
boundary is given by (|37D in both the interior and exterior coordinates. The two regions 
are connected by a common value of the boundary areal radius. 

This matching also applies to a three zone model. Suppose that there is a region 
with angular momentum for < r < fj„(t), with R = Rin and r{Rin, t) = rin{t) at the 
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outer boundary; then a vacuum (Schwarzschild-de Sitter) region for rj„ < r < font and 
for font < r an FRW region, with po = M{Rin) / Rl^^ and R = Rout at the inner boundary. 
The two boundaries evolve as geodesies in the Schwarzschild region, with m = M{Rin), 
Lin = L{Rin), El = 1 + 2E{Rin), L^ut = and = 1 + kR^^. The matching 
is achieved by following the evolution of the boundary in each set of coordinates, and 
relating common values of the boundary radius. Such a model is applicable to a universe 
with two big bangs, where the angular momentum region is born within a void in the 
expanding background some time after the main big bang. 



4. Solutions 

4.I. Solution as Elliptic Integrals for A = 

If A = 0, the polynomial inside the square root in the integrals ([TSD and (|T9D is a 
quartic. The solutions can then be written in terms of elliptic integrals, according to 
the definition of [^. The integrals are respectively: 

pr 2 

h= dr (38) 

Jro V2Er* + 2GMr^ - L^r^ + 2GML'^r ^ ^ 

r 1 

h = / dr. (39) 

Jro V2Er^ + 2GMr3 - LV^ + 2GML'^r 

The value of these integrals depends on the location of the roots of the polynomial. If 
E = 0, it can be shown that: 



Ji = — , Wr ( — — — — r + ) H sJiir; ro; 0, r_, r+) 

3V2Gm]I V 2GM ; 3(GM)f ' ' ' +^ 

'2 



Jo(r;ro;0,r_,r+) (40) 



V2GM 

Joir;ro;0,r_,r+). (41) 



V2GM 

In this, r± = L'^{1 ± a/1 - WG"^ M'^ / L"^) / AG M , and the functions J„ are given by: 

Jn{r-ro-a,b,c) = (42) 
Jro ^y\{x — a){x — bj{x — c)\ 

This result is obtained by differentiation of ^Jrir"^ — l?r I2GM + L^) with respect 
to r and integration of the result. If 7^ and the roots of r'^iE — F(r, i?)) are written 
as 0, ri, r2 and ra, the integrals become: 

_ r^E-V{r,R) 1 \GM /I 1 1 1 1 \ 
7^ +2eV 2l^ U'^'n'^'V 
/I 1 1 1 1 \ 

-^1 (43) 



2\/2E \r' ro' ri' r2' r-i 

fJo -; -; - (44) 



V2GML^ \r ro ri rg rg 
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These are derived by considering (\.{r^jE — V(r, _R))/dr, integrating the result and 
changing variable to m = in the resulting integrals. The roots rj can be obtained 
in the standard way for cubic equations (see for instance [0]). The value of the eliptic 
integral J„ depends on the number of real roots. If there are 3 roots (which may be 
repeated), the integrals are given in If there is only one real root, the term in the 
square root in J„ may be written in the form (x — a){{x — d)"^ + /^). In this case, writing 
X = a +ptan^ (0/2) with p = iy(a^^^7ip~+y^, the integrals reduce to a combination of 
elliptic integrals and elementary functions. 



4-2. Potentials with Repeated Roots 

An interesting example for non zero A is the case in which E{R) — V{r,R) has a 
triple root in r. This ensures dV/dr has a repeated root, which gives the condition 
AY = Q±{X) for X = L^/{3G^M^) and Y = 3L^. The repeated root is located at 
r = ro = 2y/XYP±{X)/3. This result and the functions Q±iX) and P±(X) are given m 
the appendix. The condition AY = Q±{X) can be rewritten as Q±{X)/X = OG^M^A. 
Given a value for A and M{R), this determines L^. For A > 0, this has solutions 
in X > 3.75 for and in 3.75 < X < 4.0 for Q_. For A < 0, this has solutions 
only for Q-{X) and X > 4.0. We shall concentrate on the former case, A > 0. In 
that case the functions Q±{X)/X have a maximum at X = 3.75, which constrains 
GM < 2/(75VA) = GM^ax- This suggests taking z = M/M^ax = z{R) as the shell 
label. X is then given by the relation Q±{X)/X = 4/(625z^), and so X = X{z) only. 
Requiring ro to be a root of E — V fixes E = V{rQ, z). The functions X, L, ro and E 
are then given in terms of A and z by the expressions 



Q±{X) 4 2 



X 625 



z^ X{z) implicitly (45) 



= = \^z^x{z) (46) 

A A 1875 ^ ' ^ ' 

ro = -^f,{z) = l^±zX{z)P^{X{z)) (47) 
V A V A 75 



^/At and r = {1 / ^/A)f {z , f) , equations ([I^ ) and ( p^D reduce to 

- I? + 2-B + + 1^ (49) 





The metric, (^), is given by ds^ = ds^/A, where 




d5' = ( 1 + ^ ) df ' - 2 ( 1 + ^ I gdfdz - r\de'' + sin' Od^P'') 
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(51) 



In this ' now denotes {d/dz)~. This solution exhibits three types of behaviour. If 
r(2;,0) = ro, the universe is static, an Einstein cluster. If the '— ' solution is taken and 
r(0,2;) < ro, the result is a universe which undergoes coasting expansion. If the '+' 
solution is used and r(0, z) > ro, the universe undergoes coasting collapse. In the latter 
two cases, the asymptotic state is a static solution. The '+' solution is only relevant to 
collapse, since E — V < for r slightly less than ro in that case, and vice versa for the 
'— ' solution. 

We may write 2{E — V{r, z)) = (f — ro)^(^ ~ fi)(f 
side of (E 



where ri : 



■Ifo ± 



^2)/(3r ) for the right hand 
^fg + 4zZV(25fi|). The solution to (||) for the 



expansion case, taking a big bang initial condition r(0, z) 



is then given by 



ri - ro 



iri 



r r 



[ro — r)(r 
273^2 



r2 



ri 



^E(5,g) 



[F(5,g) 



Denoting the integral inside (d/dz) in 

2^3 



B(f, z) 



ro[ro - r2)^/ri[ro - r2 



n(ro/(ro-r2);5, q)] . 

by B{f, z), we have also that 
[-r2n{l;6,q) + foF{6,q)]. 



(52) 



(53) 



Here 6 = arcsin ((fo - r2)f)/(fo(f - 7=2)) and q = (fo(fi - r2))/(fi(fo - f'^)). 
The functions F(0, /c), E(0, /c) and 11 (n; 0, /c) are the elliptic integrals of the first, second 
and third kinds respectively. The second result was obtained using ||2^, §3.151 4 and 
the first by writing = fof — rifo — r) and using §3.168 46 and §3.167 12. 

For the collapse case, f2 < < fi < fo < f and the two integrals become: 



r 



Bir.z 



2V3(fo - ri) 

A/fo(fi - f2) 

2v^ 



n 



^0 



ri 



^('^o - ^2) 



r2 
1 



+ 



2y3( 



A/ro(ri - r2) 



=F(i^,g) 



+ V3f^5(f , 



(54) 



f(0,^) 



=F(7,g) + 



^^2 



ro - ri 



-^E(7,g) 



V-n^2 



[^-m)(m-^) 



/r^(ro - ri) 



I'D 



— u 



(55) 



?{0,z) 



A/(l/ro - l/r2)/(l/ri - l/rs) and 

f2)(f — fo))/((fo — f2)(f — fi)). The result (E^) is obtained by 



In these, 7 = arcsin a/(u — l/fi)/(u — 1/fo) 
V = sm~^ a/ ((fi 

changing variable to -u = f^^, writing — m = (fg ^ — -u) — f^^ and using §3.131 
3 and §3.133 9. Expression ( ^4|) is derived by writing f ^ = (f — fo)^ + 2fo(f — fo) + fg 
and using ||25| §3.147 8 and §3.167 32. In both cases, the solutions have taken a common 
origin of time, so that r = for all z at t = 0. 
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Figure 3 illustrates the functions L^, fo and E for the two cases. Figure 4 illustrates 
r(f, z) for fixed values of z and f. The initial condition in the collapse case was chosen 
as r(0,z) = z^ , representing constant density. In the collapse case, expanding fo about 
z = 1 gives fo > fo(l) for 1 — z ^ 1. This is clear from figure 3c. The region in 
which dfo/dz < is not physical, as there would have to be shell crossing. Figure 4 
illustrates this, as the function r(f, z) is seen to turn over for large z at later time. This 
problem is avoided by redefining M^ax to be at the point where dfo/dz = 0, giving 
GMjnax — 2zcrit/(75VX). The point Zcrit is found numerically to be at 0.84997. 

Figure 3 also illustrates that for z 0, ro oc z^^^ in the collapse case, but ro oc 2; in 
the expansion case. Finally, we see that these universe are not closed, since E is never 
equal to —1/2. A closed model is only possible by matching onto a closed FRW exterior. 




(c) 

Figure 3 Plots of (a) LP{z), (b) E{z) and (c) fo{z) for the repeated root solution. In each plot, the solid line is 
the '+' case, relevant to collapse, and the dashed line is the '-' case, relevant for expansion. 
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Figure 4 r(f,2) at fixed z ((1), (3)) and at fixed r {(2), (4)) for the repeated root solution. (1) and (2) are 
for the collapse case and (3) and (4) are the expansion case. The values used in (1) were z = 0.001, 0.01, 0.05 to 0.5 in 
intervals of 0.05 and 0.6 to 0.9 in intervals of 0.1. For (3), these were z = 0.05 to 0.95 in intervals of 0.05. The value of z 
increases from the bottom to the top. In (2), f = 0.0 to 2.8 in intervals of 0.2, and -f increases from top to bottom. In 
(3), f = 0.01, 0.05, 0.1, 0.15, 0.4, 0.8, 1.2, 1.8, 2.5, 5.0, 10.0, 20.0 and f increases from bottom to top. 



Limiting Behaviour As r — >■ oo, these solutions approach the static Einstein cluster 
with r = ro(^). Defining D± = 2y/XP±{X)/3, X is given as = (3D^ + 5)/(4L)±), 
whence (|5])-(H) give E, L and r = ro as functions of D±. Substituting into equation 
(I), and using partial fractions, u can be evaluated as 

u= I ^/;^^ = lln(D-l) + Jln(D + l)-lnD + |ln(3D2 + l) 

— arctan v^-D + arctan | \ f—D ] . (56) 

This gives the metric of the static system in terms of the coordinates t and D. In the 
'+' case, "^0(2;) oc zs near r = and so the density e oc M'/{r'^r') is finite everywhere. 
There is no singularity in the asymptotic solution and so no naked singularity forms 
during the collapse. 

In the '-' case ro{z) oc z near z = and e oc z'"^ diverges. Using (^5])-(^) and (|56|), 
the metric near the origin is ds^ = Az^dt'^ — Bdz'^ — {^z) df2^, where A and B are 
positive constants. It is thus clear that for a radial photon, dz/dt oc z^ near the origin. 
This is integrable, so photons can escape from the central singularity - it is locally naked. 
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The singularity is always present, as a remnant of the big bang. The '-' situation is not 
appropriate to gravitational collapse, and so this naked singularity could not arise out 
of collapse. However, the presence of a naked singularity in a universe born from a big 
bang is still interesting. 

In |12[, taking M = m^R^ + m^R^ + ... and = l^R'^ + l^R^ + ... near = 0, 
they found that if Z4 7^ no central singularity forms, but one does if ^4 = 0. In this 
example X = L^/(3G'^M^), so the previous relations with I4 ^ give X oc R~^ oc M~i 
near the central singularity ai R = 0. This is the situation in the '+' case here, when 

2 

X (X z~3 ^00 near z = 0. Once again no singularity forms under these conditions. 
In the U = case, X —>■ const., which is true for the '-' case here. In fact oc z'^ 
near z = 0. The formation of the singularity in the second case and not the first is thus 
consistent with the conditions discussed by Harada et al . 



5. Summary 

A new approach to the Datta model of a spherical system of dust with angular 
momentum has been presented and investigated. Einstein's equations have been solved 
to give the metric in terms of two physical coordinates - the shell label and the proper 
time felt by the dust particles composing each shell. Conditions for regular evolution 
of these models have been considered and the possible types of evolution discussed. In 
addition to collapsing, expanding and static models, there are bouncing universes and 
oscillating universes, which may evolve without shell crossing. Some specific examples 
have been looked at in detail. The application of this approach to the formation of 
naked singularities has been discussed with reference to Einstein clusters and their 
A 7^ generalization. We find that naked singularities may form, but that these are 
not stable in the sense that spherical perturbations will lead to the singularity being 
clothed. 
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Appendix A. Classification of the Potential for non-zero A 

In the case A = 0, (p^) tells us that dV/dr = GM{r — r+)(r — r^)/r^, where 
r± = L"^ [l± yjl - 12G'2MVL2j /2GM. The potential then has 0, 1 or 2 turning 
points as X = L"^ / {SG"^ M"^) is less than, equal to or greater than 4. The signs of 
V{r±,R) indicate how many times the potential cuts the r-axis. If A is non-zero the 
potential is a quintic divided by r^. The shape of V{r, R) depends on the number of 
turning points. The transition between the various shapes occurs when the potential 
has a point of inflection, i.e. when dV/dr has a repeated 0. Differentiating equation 



(P^) and requiring dV/dr to have a double zero at r = tq, gives the condition that 
F± = Q±{X)/A - F = 0, where X = Ly{3G^M^), Y = SL^ and Q±(X) is given by: 

27(4P± - 3) 



Q±{X) 



640X2p± - 1200(1 + P±)X + 1125 



with Pi = l±^/l-ii. (A.l) 

This repeated root is at tq = XY P±{X) /3. The functions Q± are deflned only 
for X > 3.75, so that the square root in ( |A.lj ) is real. A physically reasonable solution 
must have Y c(. > 0. For A > 0, both solutions satisfy this for 3.75 < X < 4.0, and 
Q+ does for X > 4.0. If A < 0, only the Q- solution can have > 0, and X must be 
greater than 4.0. Using these functions, the shape of the potential may be classifled by 
the following scheme: 

A > : 

X < 3.75 one turning point only 

Ay < Q^{X) one turning point only 

AY = Q_{X) 2 t. p. s - inflection to left 

3.75 <X< 4.0 Q-{X) < AY < Q+{X) three turning points 

AY = Q+{X) 2 t.p.s - inflection to right 

Ay > Q^{X) one turning point only 

Ay < Q^{X) three turning points 

4.0 < X Ay = Q+(X) 2 t.p.s - inflection to right 

Ay > Q_|_(X) one turning point only 

A < : 

X < 4.0 no turning points 

Ay < (5-(X) no turning points 

X > 4.0 Ay = (5_(X) one point of inflection 

Ay > Q_{X) two turning points. 
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The cases when the point of inflection is to the left or right of the other turning 
point are distinguished by comparing tq to (3GM/(10A))^, which is where the second 
derivative of r^V{r,R) is zero. If tq > (3G'M/(10A))3, the repeated root is to the 
right of the other root, and it is to the left otherwise. This reduces to the condition 

1 2 

(10(5±)^ ~ 2X3Pj_(X)/3 < 0, which is a function of X only. For this inequality 
is satisfied for all X. For Q-, it is satisfied for X > 4.0, but not for 3.75 < X < 4.0. 

The cases that dV/dr has 3 zeros or just 1 are distinguished using the function 
{dF±/dc)a,b evaluated at F± = (in this, a = 3GM/A, b = 3L7A and c = 9GML'^/A). 
This function determines the direction the potential moves as X moves away from the 
repeated root value. In fact {—a{dF±/dc)a,b)\F±=o = Xd\iiQ±/dX + 2, a function of 
X only. This function is always negative for Q+. For it is positive ii X > 4.0 and 
negative for 3.75 < X < 4.0. 
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